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INTRODUCTION 

Extracting a broad, collimatedi ion beam from a low density 
plasma represents a difficult problem in ion-optical design. There 
have been attempts to design efficient ion extraction systems mathe- 
matically [1-5], however, the theoretical considerations involved make 

it difficult to obtain reliable solutions. The production of a broad 

.1 

ion beam is important to the development of ion thrusters designed for 
the propulsion of space vehicles. These devices obtain thrust by 
ejecting large numbers of ions at high velocities over extended periods 
of time. An accurate knowledge of the distribution of ion trajectories 
(the beam profile) is of considerable importance when applying ion 
thrusters to actual mission situations because of the thrust degrada- 
tion associated with off-axis trajectories. Broad ion beams are also 
important in a variety of ground based ibn sources. Although thrust is 
not important in these ground applications, there are other reasons for 
preferring a well collimated beam. 

This study deals with the effect of variations in extraction 
geometry and operating conditions of the extraction system on the ion 
beam obtained from an ion source. The results of this study are 
applicable to both thrusters and a variety of ground-based ion sources. 


Since this work deals only with positive ions, the term "ion" is 
used with the exclusive meaning "positive ion." 


Chapter 1 

ION BEAM FORMATION 

Grid System Operation 

The ion extraction system employed for this study was the two 
grid accelerator geometry common to most ion propulsion devices. 

Figure 1-1 depicts a portion of this multiple aperture system, showing 
the formation of an ion beamlet. The variation in electrical potential 
associated with the grid geometry of Figure 1-1 is shown in Figure 1-2. 
The grid system parameters and symbols indicated on these figures are 
listed below: 

= separation between screen and accelerator grids 
d = diameter of screen apertures 

*5 . 

diameter of accelerator apertures 
- thickness of screen grid 
t- - thickness of accelerator grid 

= Kaufman effective acceleration length [5] 
a = ion beam divergence angle 
total accelerating voltage 
= net accelerating voltage 

R = net-to-total accelerating voltage ratioT”' ~~ — 

A low density plasma is generated within a discharge chamber at a 
potential a few tens of volts above that of the screen grid, which is 
Itself usually elevated to a high positive potential with respect to 
ground [Figure 1-2], The ions produced move toward the screen grid as 
a result of random the»"mal motion and slight potential variations 
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within the plasma. At the entrance to the screen hole a stable plasma 
sheath is formed [Figure 1-1] and the tons dxit from its surface with 
the characteristic Bohm velocity [6]. 

Ion acceleration is achieved because of the potential difference 
between the plasma and the accelerator grid [Figure 1-2], Depending 
upon the amount of focusing that occurs and the subsequent ion beamlet 
diameter, dj^, most of the ions will pass through the accelerator grid 
and be expelled into the region beyond. This region is usually re- 
ferred to as the downstream side of the grids,, which for a flight 
thruster represents the vacuum of space. 

Ion Beam Divergence 

Ion thrusters, being inherently low thrust devices, must produce 
a highly collimated ion beam for efficient operation. The final tra- 
jectory imparted to each beam ion is a result of the focusing it 
receives while passing through the grids and defocus ing it experiences 
once it leaves the accelerator system. While the former mechanism ’’s 
readily understood in terms of the equipotential distribution that 
exists between an operating grid set [1-4], the factors controlling 
defocusing are less well known and warrant some discussion. 

As the accelerated ions leave the grid system, electrons are in- 
jected Into the beam to maintain a neutral charge efflux. These 
neutralization electrons, as they are called, rapidly spread through- 
out the beam producing a neutral plasma. Because of the negative po- 
tential on the accelerator grid, neutralization electrons coming 
within a finite distance of the accelerator holes are repulsed by 
electro-static forces. The result is a neutralization surface created 


slightly downstream from the accelerator grid. For multi aperture grid 
systems this surface is roughly planar [Figure 1-1]. 

Between the neutralization plane and accelerator holes the ions 
are in a region of decelerating electric field directed towards the 
negative accelerator grid [Figure 1-2]. Although the deceleration 
incurred is less than the acceleration imparted by the grid set, it is 
still sufficient to defocus the ion trajectories significantly. How- 
ever, once these ions reach the neutralization plane, electrostatic 
shielding by the beam plasma nullifies the accelerator grid's effect 
and little additional change in ion trajectories occurs. A measure of 
the overall ion beam divergence is given by the angle a [Figure 1-1]. 

It is apparent that the ion beam divergence angle is one parameter 
important in determining the ability of a particular grid geometry to . 
produce the desired ion beam colUmation. There is however another 
less obvious quantity v;hieh is of considerable interest for space 
flight applications. This is the ion beam divergence factor, f^, which 
is a ratio of the net axial thrust produced by the divergent ion beam 
to the thrust produced if the ion beam were perfectly collimated. 

These two quantities were the calculated parameters (based on experi- 
mental measurements) of primary interest to this study. 

Theoretical Considerations 

There exists a space charge flow limit to the ion current that 
can be extracted from an accelerating system of this kind. Space 
charge effects on accelerator system performance can be determined 
by solving Poisson's equation for a single aperture grid geometry 


7 


where V is the electrostatic potential and p is the charge density 
of the ions in the beam. Assuming one dimensional flow between 
parallel plane electrodes and zero electric field across the plasma 
sheath emitting surface, Equation (1-1) yields Child's current density 
law T.7] 
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Here j is the ion current density in the beam, ^ is the ion charge-to- 
mass ratio, is the total accelerating potential at a distance z 
downstream of the screen hole sheath and is the permittivity of 
free space. 

The ability of a grid set to extract the maximum ion current per 
hole for the minimum total accelerating voltage is a measure of its 
performance. This quantity is referred to as perveance and is defined 
by 
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Here the assumption has been made that the ion current J is emitted 
from a plasma sheath whose area is constant, and approximately equal 
to that of the screen hole. For this study the acceleration distance 


Z is described by the Kaufman effective acceleration length Z_ [5], 

s 

where 


d? 




(1-4) 


Actual grid set performance then can be expressed better in terms of 
a "normalized perveance per hole," which is defined as 
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Equation (1-5) indicates that a theoretical limit governs the 

I 

maximum obtainable normalized perveance per hole for a particular 
propellant. Normalized perveance per hole will be used to describe 
the ion extraction performance of each grid set investigated. 


Chapter 2 


APPARATUS 


Ion Source 

For this study, a simple mildly-divergent-field 8-cm electron- 
bombardment ion source was constructed and operated on argon propellant. 
Tungsten wire filaments were used as both the main and neutralizer 
cathode emitters. The magnetic field was derived from a long sole- 
noidal winding extending the length of the discharge chamber, with an 
additional coil winding positioned at the chamber’s rear. The geometry 
was such that the field at the front of the discharge chamber was 60% 
that of the rear. A cylindrical anode was employed and non-magnetic 
stainless steel construction used throughout. Figure 2-1 illustrates 
the basic ion source design with associated power supplies and instru- 
mentation. Further design details and operating characteristics for 
this type of source can be found in the literature [8-10], 

All source operation was conducted in a 30-cm diameter pyrex bell 
jar pumped by a 10-cm oil diffusion pump in series with a mechanical 
pump. The argon flow rate into the source was sensed by a Hastings 
mass flow meter and displayed on a digital readout. 

Grid Construction 

The screen and accelerator grids were made from thin sheet 
graphite. This material has a low ion-sputter erosion rate and a low 
coefficient of thermal expansion. In addition, because it could be 
machined easily and was readily available in the desired thicknesses, 
graphite was particularly attractive for use as a grid material. 
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The grid aperture pattern employed consisted of a nineteen hole 
hexagonal array with a center-to-center hole spacing of 2,54 mm. It 
was felt that this number of holes would adequately model the adjacent 
hole interactions found in full size grid systems. Five-centimeter 
square graphite plates containing the grid pattern were positioned on 
a masking plate which covered most of the downstream end of the dis- 
charge chamber. The array of holes covered only a small portion of the 
cross section near the axis, which insured ion extraction from a near 
uniform plasma. Calculations predicted that the greatest possible 
spatial variation in discharge chamber ion density across the grid open 
area would be less than 5 %. 

Variable grid separation was achieved by using thin mica washers 
{0.254 mm thick). These were replaced after each data run to avoid in- 
sulation breakdown and large leakage currents. The grid sets were 
fastened together using four stainless steel bolts which were properly 
insulated to prevent direct shorting. Alignment of the screen and accel 
erator grids was accomplished by hand beneath a large illuminated magni- 
■fylng glass; this straightforward technique was found to produce the 
degree of reliability desired. Grid separation could be checked when 
they were cold, but not while the source was operating and the grids 
were hot. While some grid warpage could occur during operation, the 
magnitude of this warpage should be small because the grids were carbon 
and the greatest distance between grid supports was less than 2 cm. 

Ion Beam Measurement 
(a) Beam Composition 

The ion beam produced by an electron-bombardment source is com- 
posed mainly of primary beam ions, that is ions which have received the 
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maxiinuin energy that can be imparted by the grids. However, the ion beam 
also contains a small multi-energy ion distribution produced as a result 
of primary beam ions undergoing charge exchange reactions with escaping 
neutral propellant atoms. This exchange occurs in regions between and 
downstream of the grids [11-12]. 

Charge exchange ions produced between the grids, and having suffi- 
cient energy to escape, are commonly called group 2 ions {this designa- 
tion is adopted from Kerslake et al. [11]). The acceleration they receive 
and their resultant energy depends upon where they were created between 
the grid set. Some may have energies up to hundreds of electron volts. 
Because of the improper focusing group 2 ions receive subsequent to their 
formation, most have severe off-axis trajectories and are quickly lost 
from the primary ion beam. In contrast, group 4 charge exchange ions 
originate downstream of the grids and so are created in a near field free 
region. The small electric fields which do exist within the ion beam are 
a result of plasma density variations. These fields direct the group 4 
ions further downstream, imparting to them a few electron volts of energy. 
The potential gradients are usually fairly uniform along the beam axis; 
however, they become more non-uniform and have major off-axis components 
as the beam edge is approached. The result is an ill defined beam edge 
caused by group 4 ion migration to regions outside the primary ion beam, 
(b) Faraday Probe Rake 

The basic ion beam detection apparatus used the Faraday probe to 
sense beam ion currents [13]. Typically, an individual Faraday sensor 
consists of a small flat plate electrode which is biased a few volts 
negative relative to ground and is usually positioned within the beam 
plasma at some point downstream of the ion source. Neutralization 
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electrons are reflected from the sensor's surface because of its nega- 
tive potential, leaving the ion beam impingement to be recorded as a 
net positive current. 

Distortion of the primary ion beam profile, caused by the presence 
of group 4 ions in the beam plasma, made it necessary to modify the usual 
Faraday probe design so that only primary beam ions would be detected. 

To accomplish this, a second electrode of large transparency (stainless 
steel mesh) was positioned a few millimeters in front of the Faraday 
sensors. This second electrode was held negative with respect to ground 
while the probes were biased a few volts positive above ground. 

During probe rake operation, neutralization electrons are reflected 
by the probe screen. This occurs because the mesh size, approximately 
1.2 mm square, is less than the Debye shielding distance for the maximum 
electron-ion densities found in the beam. The high energy of the pri- 
mary beam ions, several hundred electron volts, allows them to pass 
through the screen (the mesh and an open area fraction of approximately 
705i which necessitated a correction for direct interception), and impinge 
essentially undeflected on the probes. However, the much lower energy 
group 4 ions have their trajectories seriously affected by the retarding 
electric field between the screen and probe surfaces. Thus, most group 
4 ions are deflected into the screen where they are conducted away from 
the plasma. The actual effectiveness of the screened probe in suppress- 
ing group 4 ions is further discussed in Appendix I. 

In practice, twenty probes were used with a common probe screen 
electrode. The probes were mounted in line (with a center-to -center 
spacing of one centimeter) on a moveable support cradle, allowing them 
to be positioned along the beam axis. Basic construction details are 
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shown in Figure 2-2. Molybdenum was used for the ion detection surface 
of each probe because of its low secondary electron emission character- 
istics as a result of direct ion impingement. Selection of the screen 
and probe bias potentials resulted from a systematic investigation to 
determine the combination resulting in current densities which went to 
zero but did not go negative at the edges of the beam. 

Figure 2-3 illustrates the instrumentation used to monitor in- 
dividual probe currents. A digital voltmeter, measuring the voltage 
drop across a precision resistor, covered the higher ranges of probe 
current while two sensitive microammeters were used for the lower 
current ranges. Only the probe being read was biased positive, the 
surrounding probes and supporting structure were at ground potential; 
this insured minimal ion beam distortion due to secondary electron losses 
from, and group 4 ion collection on, the probe being read. 
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Chapter 3 


PROCEDURE 

Ion Source Operation 

The ion source was operated with variable beam current (controlled 
by cathode emission) at the following conditions: 


Arc voltage 

40 V 

Magnet current 

7 A 

Neutralizer coupling voltage 

0 V 

Propellant flow rate 

7 mA (equivalent) 


The arc voltage was set at 40 V to keep the production of doubly 
charged argon ions to negligible levels within the discharge chamber. 

A magnet current of about 7 ampere was found to be near optimum for the 
source used. Going beyond this value produced no additional increase 
in beam current. The neutralizer geometry consisted of a tungsten loop 
which completely surrounded the beam downstream of the accelerator grid. 
Because of the large loop size relative to the ion extraction area, effi- 
cient current neutralization could be obtained with the filament at 
ground potential. A propellant flow rate of 7 mA was the minimum which 
would give a stable discharge over the operating range of beam currents. 

Constant operating conditions (no thermal transients) was usually 
achieved within half an hour after start-up; the collection of data was 
begun soon thereafter. Beam current was controlled for the tests by ad- 
justing the refractory cathode current and hence electron mission. Each 
grid set geometry Investigated was operated over a range of emission levels 
up to the maximum beam current obtainable from the grids; the approach 
to this maximum was characterized by a rapid increase in accelerator. 


17 


18 


impingement current and a negligible increase in beam current. The 

effect of the net-to-total 

accelerating voltage ratio. 

R, on the ion 

beam profile was also investigated. Table 1 lists the 

range of voltage 

ratios used. 






Table 1: GRID 

VOLTAGES 

n » '^net 

- w 

''tot 

V 

screen 

'^accel 

V 

net 

^tot 

255 

300 

300 

600 

0.500 

375 

180 

420 

600 

0.700 

495 

60 

540 

600 

0.900 


Here v = V + V , + V 

net screen sheath arc, 


tot 


accel 


+ V 


net 


and was assumed to be 5 volts. The bell jar pressure varied 

from 5.5 X 10“^ to 8.5 x 10“^ torr during source operation. 


Obtaining a Beam Profile 

After obtaining stable source operation, the probe rake was moved 
through the ion beam with one of the central probes sensing ion current. 
The position where that probe indicated a maximum current defined the 
center of the beam. Each beam profile v/as obtained by taking a set of 
twenty ion current density measurements, corresponding to the individ- 
ual ion current received by each probe as a function of its position 
normal to the beam axis. These data were given as input to a computer 
program which calculated the ion beam divergence angle enclosing 95% 
of the total integrated beam current and the ion beam divergence factor 
associated with the entire current density profile. Further details of 
this analysis technique are discussed in Appendix II. 
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Chapter 4 

EXPERIMENTAL RESULTS 


Grid Geometry Variation 

Investigation of grid geometry variations on ion-optical perform- 

•k 

ance was done by independently varying the grid separation distance, 

accelerator aperture diameter, d^, screen grid thickness, t and 
9 as 

accelerator grid thickness, t,. To non-dimensional ize the geometrical 

a 

grid parameters they were divided by the screen hole diameter, d^, 
which remained at 0.206 cm throughout this study. Although the ion 
source was operated on argon propellant, normalized perveance is cal- 
culated for operation with mercury, since mercury is the propellant of 
principle interest for space flight thrusters. The theoretical limit 
of obtainable normalized perveance per hole, using mercury, can be 
found from Equation (1-5) and is equal to 3.03 x 10"^ amp/volt^/^. 

The reliability of the data contained within this chapter was 
verified by repeatedly testing various grid set geometries at different 
times. The results of these tests showed a maximum variation in ion 
beam divergence angle of ± 0.5 degrees and a maximum variation in ion 
beam divergence factor of ± 0.001. A tabular listing of the experi- 
mental results presented graphically in this chapter can be found in 
Appendix III. 


Ion -optical performance is defined as the degree of ion beam col- 
limation. A grid set geometry said to have "good ion-optical 
performance" produces an ion beam characterized by low beam 
divergence angles and high beam divergence factors. 
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(a) Effect of Grid Separation 

Increasing the grid separation ratio and keeping all other 
geometrical parameters constant results in significant decreases in the 
ion beam divergence angle, a, with corresponding increasts in the ion 
beam divergence factor, f^. Figure 4-1 A. Qualitatively, all curve 
trends are similar with only slight variations in the limit of obtain- 

■fr ... 

able normalized perveance per hole . In all cases the minimum beam 
divergence angle and the maximum beam divergence factor occur about a 
normalized perveance per hole of 1.0 x 10“^ amp/volt^/^. This similar- 
ity in curve shape for different values of indicates that the use of 
normalized perveance (with as the acceleration distance) is effec- 
tive in correlating performance of different geometries. Figure 4-1 B 
shows a comparison of these data with similar geometries investigated 
theoretically by Kaufman [5]. Kaufman's results agree qualitatively 
with experimental trends, but the quantitative agreement is poor. 

(b) Effect of Accelerator Aperture Diameter 

Figure 4-2 shows the effect variations in accelerator aperture 

d. 

diameter ratio ^ have on ion-optical performance. An intermediate 
grid separation ratio = 0.494) is used, with the other geometrical 

grid parameters unchanged from Figure 4-1 A. The limit of obtainable 
normalized perveance per hole increa es significantly as accelerator 


It is important to realize that constant normalized perveance does 
not mean constant current. If the total accelerating voltage 
is held constant, the beam current will increase as the grid sep- 
aration ratio is decreased [Equation 1-5]. 


ION BEAM DIVERGENCE FACTOR.fn ION BEAM DIVERGENCE ANGLE , a .(degrees) 




ION BEAM DIVERGENCE ANGLE , a , (degrees) 




R^’O.T 


THEORETICAL 

Eg/dg 

• 0.2 
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EXPERIMENTAL 

E g/ds 

O 0.247 
□ 0,494 
A 0.741 
O 1.000 
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aperture diameter ratio increases. For values of ^ above about 
0.5 the beam divergence angles and ’'iam divergence factors are inde- 
pendent of accelerator aperture size. However, for values of g— 
below 0.5 the beam divergence angles and beam divergence factors are 
dependent on accelerator aperture size. This indicates ion beamlet 

d, 

focusing becomes a function of the accelerator diameter ratio if ;q-” 

°s 

is less than about 0.5. 


(c) Effect of Screen Grid Thickness 

Figure 4-3 shows the effect variations in screen grid thickness 
tg 

ratio have on ion-optical performance. These data used an inter- 

^ .d, 

mediate accelerator aperture diameter ratio = 0.642) with the 

^s 

intermediate grid separation ratio of Figure 4-2 and the intermediate 
accelerator grid thickness ratio of Figure 4-lA. Reducing the screen 
grid thickness ratio from 0.2 to 0,1 tends to slightly increase the 
normalized perveance per hole limit, however, a further reduction has 
no significant effect on this quantity. As the screen grid thickness 
ratio decreases below 0.2, the minimum beam divergence angle and maxi- 
mum beam divergence factor are moved to slightly higher values of nor- 
malized perveance per hole. This trend is accompanied by a reduction 
in ion-optical performance in the lower ranges of normalized perveance 
per hole. Overall, there is only a minor effect on ion-optical 
performance caused by varying the screen grid thickness ratio. 

It is postulated that different screen grid thickness ratios cause 
variations in position and shape of the plasma sheath at the entrance 
to each screen aperture. These variations are thought to be the main 
cause of the small differences in the data. 
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(d) Effect of Accelerator Grid Thickness 

Figure 4-4 shows the effect variations in accelerator grid thick- 

ness ratio have on ion-optical performance. For these data, an 
^ t 

intermediate screen grid thickness ratio ( = 0.185) is used, with the 
intermediate grid separation ratio of Figure 4-2 and the intermediate 
accelerator diameter ratio of Figure 4-3. Increasing the accelerator 
grid thickness ratio increases the beam divergence angle and decreases 
both the beam divergence factor and the limit of obtainable normalized 
perveance per hole. However, these degradations are fairly uniform and 
relatively small in magnitude. 

The results shown in Figures 4-1 through 4-4 indicate that each 

geometrical grid parameter has a significantly independent effect on 

ion-optical performance. These effects can be summarized as follows. 

(i) Changes in the grid separation ratio primarily effect 

°s 

the ion beam divergence angle a and the ion beam divergence 
factor fn. 

(ii) Varying the accelerator aperture diameter ratio -A above 

°s 

about 0,5 primarily affects the limit to the obtainable 
normalized perveance per hole. 

t 

(iii) Reducing the screen grid thickness ratio ^ below about 0.2, 

Qs 

primarily moves the minimum beam divergence angle and maxi- 
mum beam divergence factor to slightly higher values of 
normalized perveance per hole. 

t 

(iv) Increasing the accelerator thickness ratio uniformly 

°s 

degrades the ion-optical performance of a particular grid 
set, but to a minor extent. 
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Variations 1n Grid Operating Conditions 

To investigate the changes in ion-optical performance for various 

accelerating voltages, selected grid set geometries were operated at 

V . 

net-to-total accelerating voltage ratios (R = vp — ) of 0.5, 0,7 and 

tot 

0.9. Figure 4-5 shows the effect of varying R on a “standard geometry" 
grid set (this designation is used because the geometrical grid para- 
meters of Figure 4-5 are similar to those currently being investigated 
on developmental ion thrusters). Increasing the net-to-total acceler- 
ating voltage ratio results in a large decrease in beam divergence 
angle and a large increase in beam divergence factor. Figures 4-6 and 
4-7 show that these changes are less pronounced for larger grid sep- 
aration ratios. As shown in Figures 4-5 through 4-7, the shape of the 
ion-optical performance curves is essentially unaltered by variations 
in t|ie parameter R. Although the curve shapes are similar, the beam 
divergence increases continuously with decreasing R. 

Figure 4-8 shows the effect of changes in R on the ion-optical 
performance of a grid set using a large accelerator aperture diameter 
ratio with all other parameters standard. The ion-optical performance 
curves are very similar to those in Figure 4-5 except that in each case 
they continue to higher values of normalized perveance per hole, an 
observation predictable from the results contained in Figure 4-2. 

Figure 4-9 shows the effect of changes in R on the ion-optical 
performance of a id set using a small screen grid thickness ratio 
with all other parameters standard. The performance curves are quali- 
tatively similar to those in Figure 4.5. However, they have reduced 
ion-optical performance at low normalized perveance per hole, with 
minimum beam divergence angles and maximum beam divergence factors 

ttanoDticKiLTrr of toe 
oate®Ai. r.' ■■ ' ' ■ 
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occurring at slightly higher values of normalized pe’^veance per hole. 
These trends are predictable from Figure 4-3. 

Figure 4-10 shows the effect of changes in R on the ion-optical 
performance of a grid set using a large accelerator grid thickness 
ratio with all other parameters standard. The ion-optical performance 
curves are similar to those in Figure 4-5 but all have been uniformly 
degraded in overall performance, trends predictable from Figure 4-4. 

In summary, varying the net-to-total accelerating voltage ratio 
produces no significant changes in the shape of the ion-optical per- 
formance curves characteristic to each geometrical grid parameter 
[Figures 4-1 A through 4-4]. Hov/ever, there is a large uniform change 
in the magnitudes of beam divergence angle and beam divergence factor 
as R is varied. For all grid set geometries investigated, reducing 
the negative potential on the accelerator grid (by going to higher 
values of R) gave no indication that neutralization electrons back- 
streamed into the accelerator system because of reduced electrostatic 
repulsion forces. 

Ion Beamlet Diameter 

The diameter of the ion beamlet between the grid sets was investi- 
gated for its dependence on normalized perveance per hole. To do this, 
accelerator grid impingement current was recorded over the operating 
range of normalized perveance for each grid set geometry studied. 
Plotting impingement current against normalized perveance per hole 
produced a curve which was essentially linear until the start of di- 
rect ion impingement where the curve rose sharply, indicating the limit 
of obtainable normalized perveance per hole. The straight line portion 


ION BEAM DIVERGENCE FACTOR, fn ION BEAM DIVERGENCE ANGLE, a, (degrees) 
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was extrapolated until the actual impingement curve departed from the 
extrapolation by an amount equal to one hundred percent of the current 
indicated by the extrapolated line (a curve deviation of fifty percent 
was also investigated, however, the results obtained were not signifi- 
cantly different from those using a one hundred percent deviation). 
Beamlet diameter was equated to the accelerator aperture diameter at 
this normalized perveance per hole condition. 

Figure 4-11 shows the effect of normalized perveance per hole on 

dh 

the ion beamlet diameter ratio for different grid separations and 

°s 

net-to-total accelerating voltage ratios. Increasing the grid separa- 

n d^ 

tion ratio ^decreases the ion beamlet diameter ratio 71 — slightly for 
dg Os 

the same value of normalized perveance per hole. Decreases in beamlet 
diameter ratio for constant normalized perveance per hole were also 
observed as the net-to-total accelerating voltage ratio was increased 
from 0.5 to 0.9. These changes were not as significant as those ob- 
tained by varying the grid separation ratio. The experimental results 
in Figure 4-11 have been compared to similar data obtained theoret- 
ically by Kaufman [5]. Kaufman’s results agree qualitatively with 
experimental trends, but the quantitative agreement is poor. 

Application to Design 

Grid sets may be designed to produce specified ion-optical per- 
formance by using the results contained in Figures 4-lA through 4-10. 
These graphical results show, that to a good approximation, complete 
independence may be assumed to exist between the effect each geomet- 
rical grid parameter has on ion-optical performance. For example, 
magnitudes of beam divergence angle and beam divergence factor can be 
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J6q 

assumed to depend solely on the grid separation ratio and the net- 
to-total accelerating voltage ratio R. Similarly, the limit of obtain- 
able normalized perveance per hole can be determined by assuming it 

d= 

depends only on the accelerator aperture diameter ratio , while the 

normalized perveance per hole at which the minimum beam divergence 

angle and maximum beam divergence factor occur can be determined (for a 

small range of normalized perveance per hole about a value of 1.0 x 

10-^ amp/volt^/^) by assuming it depends only on the screen grid thick- 
. ^s 

ness ratio -r~ . 

^s 

Figures 4-12 and 4-13 show four grid sets which were designed to 
give good ion-optical performance over different ranges of normalized 
perveance. The magnitudes of the beam divergence angles and beam 
divergence factors were predicted from Figure 4-7, the limit to the 
obtainable normalized perveance was predicted from Figure 4-2 and the 
normalized perveance at which the minimum beam divergence angle and 
maximum beam divergence factor occurred was predicted from Figure 4-3, 
Table 2 compares these predicted results with those actually measured 

'ic 

from Figures 4-12 and 4-13. Good agreement between these results is 
shown, indicating the validity of assuming complete independence be- 
tween the effect each geometrical grid parameter has on ion-optical 
performance. 

The low and middle normalized perveance designs shown in 
Figure 4-12 represent grid sets embodying the "Small Hole Accelerator 


No table of predicted performance appears for the "middle normalized 
perveance design" since these data are identical to that previously 
shown in Figure 4-7. 
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Grid" concept used to reduce the accelerator grid open area fraction 
and hence un-1onized propellant loss. This in turn improves ion source 
efficiency [14]. 


Table 2. COMPARISON OF PREDICTED AND 

EXPERIMENTAL ION-OPTICAL PERFORMANCE 


(a) Low Normalized Perveanee Design 


cl t t 

# = 1 .000, ^ = 0.469, = 0.185, = 0.370 

“s 


Nomenclature 


"^min. 

fpmax. 

N'P(fQmax.) 

N‘P 

max. 

“'"•’’max.* 


= minimum ion beam divergence angle (degrees). 

= normalized perveanee at (amp/volt^/^) . 

= maximum ion beam divergence factor. 

= normalized perveanee at f^max. 

= limit to the obtainable normalized perveanee. 
= beam divergence angle at N 
= beam divergence factor at N*P.„_„. 

lUuA 



R = 

Predicted 

Result 

0.7 

Measured 

Result 

R = 

Predicted 

Result 

0.9 

Measured 

Result 

a - 

10.1 

8.9 

6.9 

6.6 

min. 


0.63x10-^ 

0.80x10-9 

0.63x10-9 

0.80x10-9 

f Qmax . 

0.9935 

0.9943 

0.9965 

0.9963 

N*P{fpmax. ) 

0.63x10-9 

0.80x10-9 

0.63x10-9 

0.80x10-9 

N*P 

0.63x10-9 

0.80x10-9 

0.63x10-9 

0.80x10-9 

max. 


10.1 

8.9 

6.9 

6.6 

fotN-Pmax.) 

0.9935 

0.9943 

0.9965 

0.9963 
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(b) High Normalized Perveance Design 


''s 

1.000, ^ = 
°s 

t 

0.827, ^ = 

0.123, ^ = 

0.370 


R = 

Predicted 

Result 

0,7 

Measured 

Result 

R = 

Predicted 

Result 

= 0.9 

Measured 

Result 

cc . 
min. 

10.1 

10.6 

6.9 

7.5 

N*P(a . ) 

^ irnn.^ 

1.10x10-5 

1.20x10-9 

1.10x10-9 

1.20x10-9 

f Qmax . 

0.9935 

0.9931 

0.9965 

0.9960 

N*P(fpmax.) 

1.20x10-9 

1.30x10-9 

1.20x10-9 

1.20x10-9 

N-P 

max. 

2.40x10-9 

2.46x10-9 

2.62x10-9 

2.62x10-9 

“(►'•'’max.) 

13.3 

14,8 

9.8 

10.7 

foC-Pmax.) 

0.9840 

0.9864 

0.9925 

0.9916 

(c) Very High Normalized Perveance Design 



d 

1.000, ^ = 
°s 

1.000, ^ = 

0.062, = 
^s 

0.370 


R = 
Predicted 
Result 

0.7 

Measured 

Result 

R = 

Predicted 

Result 

0.9 

Measured 

Result 

“min. 

10.1 

10.6 

6.9 

7.2 

N-P'“nn-n.) 

1 .15x10-9 

1.20x10-9 

1.15x10-9 

1.30xl0~9 

f pmax . 

0.9935 

0.9931 

0.9965 

0.9963 

N*P(fpmax. ) 

1.20x10-9 

1.30x10-9 

1 .20x10-9 

1,30x10-9 

"•Pmax. 

3.30x10-9 

3.47x10-9 

3.30x10-9 

3.55x10-9 

“("•Pmax.) 

15.3 

16.8 

13.7 

13.2 

PD("-Pmax.) 

0.9780 

0.9787 

0.9865 

0.9869 


ION BEAM DIVERGENCE FACTOR, fo ION BEAM DIVERGENCE ANGLE, a , (degrees) 
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LOW AND MIDDLE NORMALIZED PERVEANCE 
GRID SET DESIGNS 


FIGURE 4-12 
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Figure 4-13 shows two grid sets designed to operate at increased 

levels of normalized perveance per hole. The high normalized perveance 

design is considered more practical of the two. It uses a smaller 

d d 

accelerator aperture diameter ratio = 0.827 compared to ^ — 1 .000) 

s ^ 

for some reduction ofun-ionized propellant loss, while employing a 
sufficiently thick screen grid to yield an acceptable ion sputter 
erosion lifetime. The very high normalized perveance grid set design 
illustrates the experimental limit to ion-optical performance. A beam 
divergence factor greater than 0.99 is evident for a normalized 
perveance of 0.4 to 3.1 x 10"^ amp/volt^/^ with a net-to-total acceler- 
ating voltage ratio of 0.9. The significant increase beyond the 
theoretical Childs' law normalized perveance limit for mercury, 

3,03 X 10"^ amp/volt^/^, indicates the approximate nature of a one- 
dimensional derivation when used for a two-dimensional configuration 
at these high operating conditions. 


Chapter 5 


A DESIGN EXAMPLE 

A Method of Predicting lon^Optical Performance 

A technique is described to predict ion-optical performance for 
developmental electron-bombardment thrusters. The technique makes ex- 
tensive use of the graphical results presented in Chapter 4. The data 
presented herein were obtained from a plasma uniformly dense across the 
ion extraction area and thus cannot be applied directly to full sized 
ion thrusters, where the plasma density within the discharge chamber is 
significantly non-uniform across the screen grid area. Incorporating 
the effects of plasma non-uniformity into calculations of ion-optical 
performance requires a knowledge of the plasma (and hence ion) density 
profile. Once this profile is obtained, the graphical data can be used 
to determine ion-optical performance by partitioning the surface of a 
full sized grid set into regions of near uniform ion density. The ion- 
optical performance is then determined for each individual region. 
Averaging these results determines the ion-optical performance for the 
entire grid set. 

To illustrate this method, the 30-cm diameter engineering Model 
electron-bombardment ion thruster, currently being developed as a 
primary propulsion source for extended space missions [15], is analyzed 
to predict its ion-optical performance. The non-dimensional ized grid 
set parameters for this thruster are as follows: 


= 0.37 


^ = 0.80 
^s 
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^ = 0.20 

''s 

= 0.27 

Developmental ion thrusters employ spherically dished grid sets to 
reduce changes in grid separation distance caused by thermal distor- 
tion [16]. However, ion beam divergence is increased when the grid 
sets are dished. This effect is compensated for by using a slightly 
different screen and accelerator aperture array [16], the ion beam 
divergence from a compensated dished grid set then closely approximates 
the ion beam divergence from a flat grid set. Hence, the results con- 
tained in Chapter 4 using a flat grid set geometry are also valid for 
compensated dished grids. 

Figure 5-1 shows the ion density profile for the 30-cm E.M. ion 
thruster. This profile was obtained from unpublished data supplied by 
Hughes Research Laboratories [17] and represents variations in ion 
density directly adjacent to the screen grid. The dependence of 
profile shape on beam current was Investigated and found to be an in- 
significant factor. The variation is therefore assumed valid for lower 
ranges of beam current. The ion density depicted in Figure 5-1 is for 
a beam current of 2.0 ampere (the designed operating level for the 
30-cm E.M. ion thruster). 

Calculations were performed for the 30-cm E.M. ion thruster to 
predict ion beam divergence factors for thruster operation at beam 
currents of 2.0, 1.5 and 1.0 ampere. To illustrate the method of 
analysis used in performing these calculations a worked example is 
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presented in step form, where the beam divergence factor is determined 
for an ion beam current of 2.0 ampere. 

Step I. 

The 30-cm grid set is partitioned into five equal annular areas 
and the average normalized ion density is determined from Figure 5-1 
for each region (calculations have shown that five equal area regions 
produce reliable results). Next, the fraction of total ion current 
that each region can extract is calculated. This is done by summing 
the normalized ion density of each region and then taking the ratio of 
the normalized ion density of each region to this sum. From these 
values, the ion current that must be extracted from each region to pro- 
duce the desired beam current of 2.0 ampere is calculated. Table 3 
shows the results obtained (region A is locateid at the grid center 
while region E corresponds to the annular area around the grid edge). 

Table 3. ION BEAM DIVERGENCE 
CALCULATIONS, STEP I 


Region 

Normalized 
Ion Density 

Fraction of Total 
Current/Region 

Current/Region 
For a 2.0 amp. Beam 

A 

1.05 

0.281 

0.562 

B 

1 .00 

0.268 

0.536 

C 

0.84 

0.223 

0.446 

D 

0.61 

0.163 

0.326 

E 

0.24 

0.065 

0.130 


Step II. 

The average normalized perveance per hole for each region is cal- 


culated by using the equation 
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N-P/H = (-^f |- 

V 3/2 d^' T 


(5-1) 


where J|^ is the ion current extracted from each region [last column 

of Jable 3] and is the total accelerating voltage (V4. = 1500 V 

2 ^ 

for the 30-cm E.H. ion thruster). The quantity (^) is found by sub- 

£ 2 

stituting the ratio -^ = 0.37 into Equation (2-4), this gives a value 

^iQ 2 ^ X 

of (^) = 0.39. The ratio j is the number of regions (X = 5) divided 

by the number of apertures in a 30-cm E.M. grid set (Y = 15173, [18]). 
With these substitutions, Equation (5-1) reduces to 


N-P/H \2.2i X 10-9 (j ) 
\ ^ 


(5-2) 


The average normalized peVveance per hole is calculated for each region 
using Equation (5-2). Table 4 lists these results. 


Table 4. ION BEAM DIVERGENCE 
CALCULATIONS, STEP II 


Region 

N-P./H (amp/volt^/^) 

A 

1.24 X 10“9 

B 

1.18 X 10"9 

C 

0.99 X 10-9 

D 

0.72 X 10-9 

E 

0.29 X 10"9 


Step III. 

The ion-optical performance expected from the 30-cm E.M. grid set 


geometry is approximated by interpolating an ion-optical performance 
curve from Figure 4-1, using the E.M. ion thruster grid separation 
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ratio = 0.37 (the E.M. ion thruster uses a net-to-total accelerating 
voltage ratio of 0.7 and operates on mercury propellant). Although the 

d- te t-a 

E.M. grid set geometry has slightly different values of ^ ^ and — 

from the values contained in Figure 4-1, these differences produce no 
significant changes in the ion beam divergence factors [Figures 4-2 
through 4-4], An ion beam divergence factor is obtained for each 
region by using the normalized perveance per hole values in Table 4 to 
read ueam divergence factors from the interpolated ion-optical perform- 
ance curve found from Figure 4-1. The results are shown in Table 5. 

Table 5. ION BEAM DIVERGENCE 

CALCULATIONS, STEP III 


Regi on 

Ion Beam Divergence 
Factor fp 

A 

0.987 

B 

0.987 

C 

0.988 

D 

0.985 

E 

0.977 


Step IV. 

The ion beam divergence factor for the entire grid set is de- 
termined by using a weighted averaging procedure. Since more current 
is extracted from the center of the grid set than the edge, the value 
of f|^ for each region is multiplied by a "weighted fraction of total 
current/region," this quantity is obtained by dividing the results in 
column three of Table 3 by the fraction of total current/region in 
region E. These results are then averaged in the usual manner to 
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u>i>- J' 


give the total ion beam divergence factor. Table 6 illustrates the 
procedure , 

Table 6. ION BEAM DIVERGENCE 
CALCULATIONS, STEP IV 


Region 

Weighted Fraction of 
Total Current/Region 

Weighted Beam 
Divergence Factor 

Average 
f f 

A 

4.3P3 

4.267 


B 

4.123 

4.069 

^ 15.172 

■^D 15.384 

C 

3.430 

3.389 


D 

2.508 

2.470 

fp= 0.986 

E 

1 .000 

0.977 


Total = 15.384 

Total = 15.172 



Similar calculations were also performed for an ion beam current 
of 1.5 and 1.0 ampere. Table 7 lists the results of these calculations. 


Table 7. ION BEAM DIVERGENCE CALCULATIONS, 
1 .5 AND 1 .0 AMPERE ION BEAM 


Region 

N-P/H (amp/voltV^) 

fo 

Weighted fp 

Average fp 






A 

0.93 X 10-5 

0.988 

4.271 

^ _ 15.160 
^D 15.384 

B 

0.89 X 10-5 

0.987 

4.069 

C 

0.74 X 10“5 

0.986 

3.382 


D 

0.54 X 10-5 

0.982 

2.463 

fp = 0.985 

E 

0.22 X 10“5 

0.975 

0.975 


^beam"^ 





A 

0.62 X 10-5 

0.984 

4,254 

. _ 15.100 
15.384 

B 

0.59 X 10-5 

0.983 

4.053 

C 

0.50 X 10-5 

0.982 

3.368 


D 

0.36 X 10-5 

0.978 

2.453 

fp = 0.981 

E 

0.15 X 10-5 

0.973 

0,972 



io i^UOE 


rnr 

V 'i; . : 


*1 » <» 
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Figure 5-2 shows a plot of ion beam divergence factor versus ion 
beam current for the 30-cm E.M. ion thruster. The predicted results 
from Tables 6 and 7 lie within the data spread obtained from independ- 
ent experimental determinations. The calculated results are considered 
most reliable because they are based on data where variables such as 
grid separation were more easily controlled. 

It is interesting to note that the grid set used with the 30 -cm E.M, 
ion thruster is operating below its normalized perveance per holft limit. 
Figure 4-2 indicates an obtainable normalized perveance per hole limit 
of 2.3 X 10“^ amp/volt^/^ for an accelerator aperture diameter ratio of 
0.8. At a beam current of 2.0 ampere the average normalized perveance 
per hole is 0.89 x 10"^ amp/volt^'^^, which is also somewhat below the 
Optimum value for ion-optical performance of 1.05 x 10" [Figure 4-3]. 

For this beam current the highest current grid apertures are operating 
at a normalized perveance per hole of only 1.29 x 10"^ amp/voit^/^. 


Reducing Accelerator Grid Effective Open Area 

Grid sets presently being used in developmental ion thrusters have 
Small Hole Accelerator Grid (S.H.A.G.) optics. As was previously men- 
tioned in Chapter 4, reducing the accelerator aperture diameter ratio 

TP reduces the loss of un-ionized propellant and improves ion thruster 
°s 

efficiency. However, Figure 5-3 illustrates that reductions in the 

d, 

ratio ;p significantly decrease the fraction of obtainable Childs' law 

^ dg 

current density. Figure 5-3 was obtained by varying ^ for a constant 

s, ^ 

grid separation ratio ^ , the resultant changes to the limit of ob- 
tainable normalized perveance can then be thought of as changes in the 


obtainable Childs' law current density for that grid separation ratio. 




ACCELERATOR APERTURE DIAMETER RATIO, dg/d 
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£g/ds 

o 0.494 

□ 1.000 
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FRACTION OF OBTAINABLE CHILD'S LAW CURRENT DENSITY 


COMPARISON OF PREDICTED AND EXPERIMENTAL 
ION BEAM DIVERGENCE (30-cm E.M. ION THRUSTER) 

FIGURE 5-3 


It was assumed the theoretical normalized perveance per hole limit for 
mercury, 3.03 x 10"^ amp/volt^^^, corresponded to an accelerator 
aperture diameter ratio of 1.00 (data were obtained which exceeded this 
theoretical limit [Figure 4-13], however, the discrepancy is thought to 
be a result of the approximate nature of the one dimensional Childs' 
current density derivation when used for a two dimensional con- 
figuration). 

Although accelerator aperture diameter does not appear in Childs' 
current density law explicity [Equation 1-2], Figure 5-3 reflects the 
fact that reducing the ratio ^ causes direct ion impingement on the 
accelerator grid at lower ion extraction current densities. For ex- 
ample, decreasing the accelerator aperture diameter ratio from 1.000 
to 0.642, using the same grid separation ratio, reduces the obtainable 
Childs' law current density by half [Figure 5-3]. This substantial re- 
duction in ion current density is a detrimental factor in the use of 
small hole accelerator grids. 

The un-ionized propellant atoms which escape through the acceler- 
ator grid are in free molecular flow. Under this flow condition, the 
effective open area each accelerator aperture presents to the escaping 
propellant atoms is less than its physical open area. The extent of 

this effective area reduction depends upon the ratio of accelerator 

t 

grid thickness to accelerator aperture diameter, , and is quantita- 

°a t 

tively defined by the Clausing factor [22]. Increasing the ratio 
decreases the accelerator grid effective open area. Furthermore, 
previous results show that increasing the accelerator grid thickness 
(for a constant accelerator aperture diameter) only slightly decreases 
the limit to the obtainable normalized perveance per hole [Figure 4-4], 


EEPBODUCIBTU'PI OF THM 
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As a result, if a set reduction in effective open area is desired, the 
ion current density which can be extracted using a grid with slightly 
decreased accelerator aperture diameter and increased accelerator grid 
thickness is greater than the ion current density which can be ex- 
tracted using a grid with reduced accelerator aperture diameter only. 
Table 8 compares this alternate method for reducing the effective ac- 

s 

celerator grid open area with the S.H.A.6. optics presently being 

investigated for the 30-cm ion thruster [23] {the effective reduction 

in open area and ion current density is determined relative to the 

d 

standard E.M. accelerator grid geometry whose parameters are = 0.80 
and ^ = 0.27). 


Table 8. ALTERNATE METHOD TO REDUCE THE ACCELERATOR 
GRID EFFECTIVE OPEN AREA 

Parameter Effective Reduction Ion Current Density 

Varied from In Open Area from Reduction from 

Standard E.M. Standard E.M. Standard E.M. 

Accelerator Grid Accelerator Grid Accelerator Grid 


S.H.A.G. ^-.0.80^0.60 48% 41% [Figure 5-3] 

Optics S 

d 

Alternate 0.80 ^ 0,70 47% 33% [Figures 5-3 

Method '^s and 4-4] 

ta 

0.27 0.67 


Table 8 shows that 12% more ion current can be obtained by using 
the alternate method instead of the S.H.A.G. optics to achieve a set 
reduction in effective accelerator grid open area. 

Another advantage to using a thicker accelerator grid is longer 
accelerator grid lifetime because more grid material is available to 
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resist ion-sputter erosion [24], The S.H.A.G, optics have shown an 
electron backstreaming limit occurring at higher net-to-total acceler- 
ating voltage ratios than obtainable with the standard E.M. accelerator 
grid geometry [14]. An increased electron backstreaming limit would 
be expected also for a thick accelerator grid, where the negative ac- 
cellerator potential is more uniformly distributed across each 
aperture. 
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CONCLUSION 

A detailed experimental investigation has shown that to a good 

approximation, each geometrical grid parameter independently affects 

one aspect of ion-optical performance. These trends are listed below. 

(i) The magnitudes of ion beam divergence angle and ion beam 

divergence factor are controlled by variations in the grid 
A 

separation ratio (at a constant net-to-total accelerating 

voltage ratio and normalized perveance per hole). 

(ii) The limit to the obtainable normalized perveance per hole is 

controlled by variations in the accelerator aperture diameter 
d 

ratio . 

•^s 

(iii) The normalized perveance per hole at which the minimum beam 

divergence angle and maximum beam divergence factor occur is 

controlled {for a small range of normalized perveance per hole 

about a value of 1.0 x 10"^ amp/volt^^^) by variations in the 

tg 

screen grid thickness ratio ^ . 

s t 

(iv) Variations in the accelerator grid thickness ratio, , 
slightly change the overall ion-optical performance. 

A graphical technique has been developed incorporating these 
results which can be used to predict ion-optical performance for an 
arbitrary ion source and grid set geometry combination. 

The graphical results are quite general and may be used for any 
application where a two-grid optics system is used to produce a broad 
ion beam. 
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APPENDIX I 

Figure A-1 illustrates the effectiveness of using a screened probe 
rake for removing low energy charge exchange ions from the primary ion 
beam. The measured ion beam divergence angle, ot, increases by an average 
of 25% and ion beam divergence factor, fp, decreases by an average of 
0.3% when using an unscreened probe rake. These data were obtained for 
the high normalized perveance design grid set [Figure 4-lA] but are 
representative of the trends expected with other grid set geometries. 

The most significant difference in the beam profile when using an un- 
screened probe rake was the absence of a zero in ion current for the 
outer probes. This indicates the large extent of the background charge 
exchange plasma surrounding the primary ion beam. 


EFFECT OF SCREENED FARADAY PROBES ON THE ION BEAM PROFILE 

FIGURE A-1 


l-Or UNSCREENED PROBES 



PROBE NUMBER 


NORMALIZED PERVEANCE = 0.38 X 1 0“® ( amp/volt^-^) 

UNSCREENED SCREENED 

a 15.9° 13.5° 

fu 0.9852 0,9880 


NORMALIZED PERVEANCE = I .OOX IO"^(amp/volt^/ 2 ) 

UNSCREENED 

a 13.5° 
fj3 0,9901 


SCREENED 

10.3° 

0.9917 

Ql 


NORMALIZED PERVEANCE = 1 .50X IO“®{amp/volt^/ 2 ) 

UNSCREENED SCREENED 

a 14.2° 10.9° 

fo 0.9893 0.9930 
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APPENDIX II 


The ion beam divergence angle, a, and the ion beam divergence 
factor, fp, were calculated from the probe rake data assuming an ion 
beam point source and circular symmetry for each ion beam profile. The 
following equations were used in their determination 


i\ 

a = Arctan 


77 


I J . ,R. Cos0. 
_ i^l ^ ^ ^ 


77 


I J,- Rt ^ 

i=l ^ ^ 


where 

= radius normal to the beam axis defining a cone enclosing 

95% of the total integrated beam current {the probe rake was 

positioned 17.5 cm downstream of the ion source) 

J. = probe current for the i^*^ probe (an interpolation routine 

was used to increase the effective number of probe data 

points from 20 to 77) 

th 

R. - distance from the i probe to the ion beam axis 

0 - = angle formed by intersection of a line extending from the 

th 

center of the ion source grid set to the i probe and the 
ion beam axis. 




?GOa 
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APPENDIX III 


Experimental Results 
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0.875 

16.47 

: -.21 
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16.19 
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5.353 

1 .665 

16.95 
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20.37 
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1.599 

0.790 

19.03 

1.887 

0.932 

18.75 

2.190 

1.082 

18.38 

2.494 

1.232 

18.13 


Both perveance and normalized perveance per hole have been 
for use with mercury propellant. 


0.9682 
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0.9849 

0.9843 

0.9826 

0.9793 

0.9686 

0.9724 

0.9748 
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0.9801 

corrected 
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